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ABSTRACT 

The freshwater genus Biomphalaria Preston, 1910 is a well 
studied genus of freshwater snails, due to its importance as an 
intermediate host of the human parasite Schistosoma mansoni 
(Platyhelminthes). Soft parts are used to differentiate among 
modern species of Biomphalaria , but these parts are not pres¬ 
ent in fossils. The aim of this paper is to show another approach 
to discriminate among Biompahalaria species, by means of the 
analysis of the logarithmic spiral of the shell. We compared 
five species, two modern and well known ones, Biomphalaria 
peregrina (d’Orbigny, 1835) and Biomphalaria tenagophila 
(d’Orbigny, 1835), a fossil one, Biomphalaria walteri (Parodiz, 
1969), and two new fossil species. All fossils are from the 
Queguay Formation (Late Cretaeeous, Uruguay). The vari¬ 
ance of the spiral was analized using a Kruskal-Waliis test with 
a posteriori Mann-Whitney pairwise comparisons. A morpho¬ 
logical analysis had to be done to determine whether the 
specimens belong to different species. Our conclusions are 
that the analysis of the variance of the spiral coefficient can 
be used as a complementary character, but not as a primary 
method of differentiation between species. A new species is 
described as a result of the combined results of the statistical 
analysis and the traditional morphological description. 

Additumal Keywords: Late Cretaceous, logarithmic spiral, 
freshwater fossils 


INTRODUCTION 

Species of the freshwater gastropod Biomphalaria Preston, 
1910 are distributed in the tropics and subtropics of the 
Americas and Africa (Taylor, 1988). Several species have 
medical importance because they are intermediate hosts of 
the human parasite Schistosoma mansoni (Platyhelminthes), 
whieh affects more than 200 million people in the world 
(WHO, 1985; Crompton et ah, 1999; Cliitsulo et ah, 2000). 

The shell in the genus is rather simple and not very 
informative taxonomically. Shell morphology of the 


genus, therefore, is used only as a secondary taxonomic 
criterion. Taxonomy of the genus is based on the anatomy, 
mainly the reproductive system (Paraense and Deslandes, 
1959; Paraense, 1975, 1981, 1984, 1988; Luz et ah, 1998), 
or even, and most recently, on molecular genetics, in par¬ 
ticular DNA (Vidigal et ah, 2000, 2004; Spatz et ah, 
1999; Dejong et ah, 2001). A relatively recent review 
of the systematics of the genus can be found in Janie 
et al. (2011). 

Some authors did statistical analyses of some shell 
dimensions (height, width), or number and shape of the 
whorls, but this methodology gave only approximate 
results (Jarne et ah, 2011). Identification of species based 
only on the shell morphology is a common problem in 
mollusks, in particular in gastropods, and non-traditional 
morphometric solutions has been tried. For example, 
to distinguish between species of Physa , Samadi et al. 
(2000) studied shells and anatomical characters of some 
populations. They estimated Raups (1966) parameters 
and analysed the variance. Unfortunately, the land¬ 
marks required in this method cannot be obtained for 
a planispiral shell, consequently the procedure is not 
useful for Biomphalaria. Johnston et al. (1991) modeled 
the spiral growth using geometric morphometries, but 
they used shells of the marine genus Epitonium that has 
a conispiral shell with axial varices along the shell, and 
they could use the varix-suture intersections as land¬ 
marks. Determination of landmarks requires homologous 
intersections of lines, necessarily constant and repeatable 
in all specimens (Bookstein, 1991; Park et al. 2013). Due 
to its simple, planispiral shell, Biomphalaria does not have 
such intersections. 

Therefore, neither Raups parameters nor geomet¬ 
ric morphometries analyses can be done in shells of 
Biomphalaria. In this paper we tried a new approach for 
distinguishing speeies with planispiral shell. The shape of 
the spiral is a primary feature being used today to distin¬ 
guish a priori different species of Biomphalaria (Bonetto 
et ah, 1982, Johnston et al., 1991). In the present work, 
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we quantify the spiral curve, in order to test its useful¬ 
ness in taxonomic differentiation among several fossil 
and modern species. 

Tin *ee fossil and two modern species were considered 
(Figures 1-10). Fossil shells of Biomphalaria walteri 
(Parodiz, 1969), and two fossil species assigned to 
Biomphalaria , one of them described at the conclusion of 
the paper as a new species, come from the limestones of 
Queguay Formation (Dite Cretaceous) in Uruguay. One 
of tlie undetermined fossils was previously unknown (our 
Biomphalaria sp. 1), and another erroneously allocated to 
Scolodonta semperi Doering, 1874 by Morton and Herbst 
(1993). This latter species is herein described as new. The 
modem species used for die analysis were B. peregrina 
(d’Orbigny, 1835) and B. tenagophila (d’Orbigny, 1835). 
They currently inhabit Uruguay and surrounding regions. 


MATERIALS AND METHODS 

The diree fossil specimens mentioned above, Biomphalaria 
walteri (Figures 1, 2), Biomphalaria sp. 1 (Figures 3, 4), 
and Biomphalaria reverse new species (Figures 5, 6) 
were compared with two recent species, Biomphalaria 
peregrine (Figures 7, 8) and Biomphalaria tenagophila 
(Figures 9, 10). The fossil specimens of Biomphalaria 
walteri and part of the specimens of Biomphalaria sp. 1 
belong to the paleontological collection of Facultad 
de Ciencias (Montevideo, Uruguay) FCDP (4668, 6462, 
6457, 6463, 6465); the rest of the specimens of 
Biomphalaria sp. 1 (no collection number) and the speci¬ 
mens of Biomphalaria reverse new species (PZ--CTES 
5345) belong to the Paleontological Collection of Uni- 
versidad Nacional del Nordeste (Corrientes, Argentina). 



Figures 1-10. Shells of Biomphalaria species. 1-2. Biomphalaria walteri (Parodiz, 1969); 3 4. Biomphalaria sp. 1. 5-6. Biomphalaria 
reversa new species, holotype. 7-8. Biomphalaria peregrina (d’Orbigny, 1835). 9-10. Biomphalaria tenagophila (d’Orbigny, 
1835). For all species, the first figure is in apical view, second figure umbilical view. Seale bar=5 mm. 







F. Cabrera et al., 2016 


Page 7 


The recent and Quaternary specimens of B. peregrina and 
B. tenagophila are hosted in the collections of Museo 
Nacional de Historia Natural (MNHN, Montevideo, 
Uruguay, no collection number); Universidad Federal 
de Santa Maria (Santa Maria, Brazil) (BR100888, 101084, 
101888, 100747); Paleontological Collection of Facultad 
de Ciencias, (FCDP 2097), and in a private collection 
(Gustavo Lecuona collection, GLC, without collection 
number). In total, 143 specimens were used in this work: 
30 of B. walteri; 14 of Biomphalaria 1; 33 of Biomphalaria 
reversa new species; 36 of B. tenagophila, and 34 of 
B. peregrina. 

The 143 specimens were photographed, digitized, and 
the spirals highlighted. Rectangles corresponding to each 
of the whorls of the spiral were plotted. Next, using 
the program Image Tool v. 3.00 (Wilcox et al., 2002), we 
measured height and width of the last three whorls of 
each of the shells from outside in (Figure 11). We used 
only the last three whorls because in most fossil specimens 
the initial whorls were broken or had encrustations. 

The value of the logarithmic spiral coefficient (0) was 
calculated as: 0 ^ AB/AD , AB being the width of a whorl, 
and AD the height; although this may not the exact value 
of 0, it may be a good approximation. The values of 0, 0', 
and 0" were obtained for each specimen, with AB/AD , 
A'B'/A'D and A"B ,, /A"D" being the rectangles of the 
last three whorls (Figure 11). To obtain a unique value 
of 0 per specimen the values were centered and stan¬ 
dardized. Next, absolute 0 values were converted to their 
natural logarithm to make them more comparable. 

The selection of parametric or nonparametric analy¬ 
sis of variance (ANOVA) for analyzing the values of 0 
between species was determined via Shapiro-Wilk tests 



Figure 11. Step by step process to obtain 0 values of the 
logarithmic spiral of the planispiral shell. 0 can be approached 
by the ecuation 0 ~ AB/AD. The 0 values of the three lasts 
whorls were averaged to obtain an unique value of 0 for each 
specimen. Species: Biomphalaria peregrina. Scale bar=5 mm. 


(SW) for normality of data and by Levenes tests (LV) for 
homogeneity of the variance. When deviations from nor¬ 
mality and homogeneity were detected, a non-parametrie 
Kmskal-Wallis test with a posteriori Mann-Whitney pair¬ 
wise comparisons (Bonferroni corrected) test were 
applied (Sokal and Rohlf, 1998) to compare 0 values of 
each species. 

Then, height, width, and 0 of each specimen were ana¬ 
lyzed using Multivariate Analysis of Variance (MANOVA) 
and Canonical Variate Analysis (CVA) (Manly, 1994). The 
CVA produces a scatter-plot along the first canonical axis, 
showing the maximum separation between groups (multi¬ 
group discriminant analysis). The axes are a linear combi¬ 
nation of the original values, and the eigenvalues indicate 
the variation of the axis and which variables contrib¬ 
ute the most to discriminate the groups. The discriminat¬ 
ing power of the variables was evaluated using Pillai 
trace because this may be more robust than Wilks lambda 
(Rencher, 2002). The relationship between width and 
height was analyzed using a linear regression model (Sokal 
and Rohlf, 1998), and the hypothesis of isometric growth 
(b = 1) was assessed applying Student t-test. PAST 
(v. 2.17c) (Hammer et al., 2001) was used for all statistical 
analysis. For all tests, the significance level used was 
p= 0.05. 

RESULTS 

The regression model (Figure 12) between width and 
height was positive and significant (r = 0.99867, n=143, 



Height 

References 

+ Biomphalaria walteri Biomphalaria peregrina 

m Biomphalaria 1 X Biomphalaria tenagophila 

O Biomphalaria reversa sp. nov 

Figure 12, Regression model between width and height; 
r=0.99867, n=143, p<0.05, b=1.0008, s b =0.0043523. 











Page 8 


THE NAUTILUS, Vol. 130, No. 1 


B. peregrine 



B tenegophlle 



B. reverse 
nov. sp. 


h 


Biomphalaria 1 


B. welieri 


■ 0 - 

HD 



0 0.8 1.6 2.4 3.2 4 4.8 5.8 6,4 

Thela 


Figure 13* Box-plot of spiral coefficients (0 values) for each 
species. The box represents 25 and 75 percent quartiles, 
respectively. The median is shown with a horizontal line inside 
the box. The minimal and maximal values are shown with short 
horizontal lines (“whiskers”). 


p < 0.05) and showed an isometric shell growth when 
all species were considered (Figure 12) (b=1.0008, 
Sb^O.0043523 and p=0.85528); an allometric effect is 
therefore discarded. 

A boxplot was performed to depict the basic statistical 
values of 0 (Figure 13). The extant species Biompholaria 
peregrina shows a greater dispersion than B. tenagophila , 
and this last one is bound to a restricted area of the 
graphic. Regarding the fossil species, B. walteri and 
Biompholaria sp. 1 have almost the same dispersion, while 


Biomphalaria reversa new species is dispersed to the right 
of the chart and only shares part of the minimal values 
with the maximal values of B. walteri. 

The 0 data were non-normal (SW=0.94; p=9.64e-06), 
non-homogeneous among species (LV=9.8e-18; p=2.Xe^38), 
and varied significantly among the analyzed species 
(Kru shall-Wallis, KW=98.52, p=2.Q34e-20). According to 
the Mann-Whitney (MW) a posteriori test with Bonferroni s 
correction, the modern species B. tenagophila and 
B. peregrina show differences between each other 
(psame= 1.879E-07). There is a similarity between the 
fossils B. walteri and Biomphalaria sp. 1 (p sa me = 0-053)> 
and significant differences between B. walteri and 
Biomphalaria reversa new species (p same =4.326E-10) 
and between Biomphalaria sp. 1 and B. reversa new 
species (p same =8.234E-07). Biomphalaria tenagophila 
differs from all fossils as follows: B. walteri (p same = 
L265E-07), Biomphalaria sp. 1 (p same =0.00349), and 
B. reversa new species (p same =1.003E-ll). However, the 
fossils B. walteri and Biomphalaria sp. 1 are both similar 
to B. peregrina (p S ame=l with B. walteri , and p sa me= 
0.5528 with Biomphalaria sp. 1) (Table 1). 

MAN OVA Hotelling’s p values (Bonferroni corrected) 
(Table 2) show significant differences between Biomphalaria 
tenagophila and B. peregrina (p=5.50382E-10); B. walteri 
shows differences with the remainder of the study species 
(p=0.0491237 with Biomphalaria sp. 1; p=2.35686E-16 
with Biomphalaria reversa new species; p=4.02082E-08 
with B. tenagophila; and p=0.000277502 with B. peregrina); 
Biomphalaria 1 shows differences with the other fossils 
but not with the extant species (p=2.69672E-12 with 
Biomphalaria. reversa new species; p=0.0522732 with 
B. tenagophila , and p=Q. 107649 with B. peregrina); 
Biomphalaria reversa new species shows differences 
with all die species (p=i.50006E-24 with B. tenagophila 
and p=3.53252E-13 with B. peregrina). CVA was signifi¬ 
cant (Filial trace=0.9918; F=17.53, p<<0.05). The first 
factor explains 89.55% of the total variability, and the 


Table 1. Kruskal - Wallis test for 0 values H (Chi 2 ) = 98.52; He (tie corrected) — 98.52; p same — 2.034E-20. Mann-Whitney 
pairwise comparisons, Bonferroni corrected (below diagonal) \ uncorrected (above diagonal). 
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Table 2. MANOVA/CVA pairwise comparisons for width, height and 0. Hotelling’s p values, uncorrected significance. It shows 
significant differences between most specimens (p-values<<a 0.05). 
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second the 6.2%. The variables that most contributed to 
distinguish among species groups were: height in the first 
factor and width in the second factor. CVA scatter-plot 
shows a well-separated Biomphalaria reversa new species 
in the left of the graphic, and Biomphalaria sp. 1 and 
B. tenagophila in the right. Biomphalaria waited has a 
wide distribution in the right zone, and B. peregnna is 
distributed throughout the chart (Figure 14). 

Analytical Discussion: The MW and MANOVA anal¬ 
yses for Biomphalaria peregnna and B. tenagophila 
show significant differences between the two species 
(Tables 1, 2). The CVA scatter-plot shows B. peregnna 
as having a wide distribution on the chart. This species 
shows large intraspecific variability, also recognized at 
the molecular level (Caldeira et al., 2001; Carvalho 
et al., 2001), and overlaps with tire other extant species, 
B. tenagophila . Significant differences for 0 and for the 
other morphological variables suggest that with a simple 
shape as the logarithmic spiral of Biomphalaria , analyses 
of specimens of distinct species can yield similar shell 
shape results. 

The statistical analysis of the extant species allows 
us to compare the statistical behavior of the fossils. 
Biomphalaria waited and Biomphalaria sp. 1 do not 
show statistical differences, the MANOVA p-value for 
both is almost 0.05 (Table 2), and both show almost in 
the same area of the scatter plot, so there are no reasons 
to consider Biomphalada sp. 1 as a different species from 
B. waited . Regarding B. reversa new species, this fos¬ 


sil species presents significant differences with both 
B. waited and Biomphalada sp. 1. Moreover, in the scat¬ 
ter plot, the new fossil species presents a well-separated 
distribution in respect to the other two fossil species. 

When the two extant species are compared with the 
fossils, 0 value for B. peregnna is similar to the values 
for B. waited and Biomphalada sp. 1 and MANOVA 
test also showed similitudes to Biomphalada 1 with the 
two extant species. Biomphalada reversa new species 
showed significant differences with the two modern spe¬ 
cies in both statistical analysis. In the CVA scatter-plot, 
B. waited and Biomphalada 1 overlap with the two 
extant species, and Biomphalada reversa new species 
overlaps with part of the distribution of B. peregnna. 

The pairwise comparisons (MW) and the multivariate 
analysis of the variance (MANOVA) show that the spiral 
growths of B. waited and Biomphalada sp. 1 are almost 
coincident, and both differ from that of Biomphalada 
reversa new species. When observing the results for dis¬ 
tribution of multigroup discriminant analysis for the three 
fossils species in the CVA scatter plot, Biomphalada 
reversa new species is well separated from the other two, 
and these overlap with each other. When the entire spe¬ 
cies set is observed in the CVA scatter plot, all of them can 
be considered as B. peregnna , so the scatter-plot by itself 
is not conclusive. We need to consider that B. peregnna is 
a species that shows a well-known, broad intraspecific 
variation, and as we used samples from Brazil and from 
Uruguay, perhaps regional variation is the reason for its 
wide distribution in the chart. Also, as mentioned before, 



■ Biomphalaria 1 

Axis 1 

O Biomphalaria reversa sp. nov 
C Biomphalaria peregrin a 
X Biomphalaria tenagophila 

Figure 14. Canonical Variate Analysis (CVA) scatter plot for height, width and 0 values for all the species. 
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it is expected that a simple shape as it is a logarithmic 
spiral could be repeated in the growth shape of differ¬ 
ent species. 

To summarize, the statistical results show that the 
three variables explain the growth of each group as a 
whole. When it was expected that only 0 would explain 
the variation, it turns out that the three variables 
together explain it. Hence, 0 is a variable that is involved 
in the relative shell growth of each species, and not a 
variable that could explain the growth shape by itself. It 
is for that reason that the statistical discrimination may 
be complemented with qualitative moi*phological char¬ 
acteristics. This evidence leads to the formal description 
of a new species the new species as follows. 

SYSTEMATICS 

Superfamily Planorboidea 
Family Planorbidae 
Subfamily Planorbinae 

Genus Biomphalaria Preston, 1910 

Type Species: Biomphalaria srnithi Preston, 1910 by 
monotypy. 

Biomphalaria reversa new species 

(Figures 5, 6, 15-17) 

Scolodonta semperi. —Morton and Herbst, 1993, p. 450, 
pi. I, figs. 9-10) (non Doering, 1874) 

Diagnosis: Shell very small, sub-circular in apical 
view, oblong in apertural view; dextral spire with four 
whorls. In apical view, each whorl covers part of the 
previous whorl. Two last whorls visible in abapical view. 

Description: Shell very small, sub-circular in apical 
view, planispiral, oblong in apertural view; dextral spire 
with four convex whorls. In apical view, early whorls 
closely stretched, being looser and wider in last portion 
of spire; last whorl is larger than preceding ones, cover¬ 
ing partially previous whorl. In umbilical view, only last 
two whorls visible. Umbilicus present. 

Type Material: Holotype: FCDPI 7230 (Figures 5, 6, 
15-17), width 6.81, height 5.79, from type locality; 
Paratypes, three specimens, FCDPI 7830, segregated 
as 7830a (Figures 18-20); 7830b (Figures 21-23); and 
7830c (Figures 24-26); all from Piedras Coloradas, 
Paysandu Department, Uruguay (Quegay Formation, 
Late Cretaceous) 

Type Locality 7 : Palmar, Soriano, Paysandu Depart¬ 
ment, Uruguay (Quegay Formation, Late Cretaceous). 

Other Material Examined: FCDPI 7254, 7267, 7271 
(4 specimens) Trinidad, Flores Department; 7220, 7253, 
7259 (107 specimens) Piedras Coloradas, Paysandu 
Department; 7241, 7261, 7265, 7266, 7276, 7277 (78 
specimens) Quebracho, Paysandu Department; 7038, 


7039 (141 specimens) Nearby Algeria Town, Rio Negro 
Department; 7829 (65 specimens) Palmar, Soriano 
Department; PZ-CTES 5345 (53 specimens) Piedras 
Coloradas, Paysandu Department (All Queguay Forma¬ 
tion, Late Cretaceous). 

Geographic and Stratigraphic Ranges: Nearby 
Algorta Town, Rio Negro Department; Trinidad, Flores; 
Palmar, Soriano; Quebracho Town and proximities of 
Piedras Coloradas Town, Paysandu Department, Uruguay 
(All Queguay Formation, Late Cretaceous). 

Etymology: Latin reversa , reverse, in referenee to the 
dextral shell orientation, in opposition to the sinistral 
orientation of the closest fossil species. 

Comparative Remarks: The specimens allocated to 
Scolodonta semperi by Morton and Herbst (1993) repre¬ 
sent Biomphalaria reversa new species. These speci¬ 
mens do not fit the description of Scolodonta semperi 
nor of any Scolodonta for that matter (Doring, 1875; 
Hausdorf, 2006). Species of Scolodonta have a discoidal 
shell with slowly increasing whorls and slightly raised 
spire (together with soft parts characters that are obvi¬ 
ously only observable in recent species.) None of these 
characteristics are present in the fossil specimens, their 
characters closely resembling those of Biomphalaria 
(see below). 

Biomphalaria reversa new species has its spiral growth 
very similar to B. ivalteri. However, it is clearly smaller 
than B. ivalteri. Additionally, B. reversa new species is 
dextral, and B. ivalteri is sinistral. Biomphalaria reversa 
has less whorls than B. ivalteri as well. In apertural view, 
B. reversa new species is more oblong and flattened than 
B. ivalteri. 

The recent species Biomphalaria peregrina is similar 
to B. reversa new species in spiral growth and general 
shell shape (Figures 15-32). Additionally, both are dex¬ 
tral. However, B. peregrina is larger than B. reversa new 
species, and B. peregrina has more than six circular 
whorls with well-defined sutures, while B. reversa has 
four slightly convex whorls, with barely visible sutures. 

Biomphalaria tenagophila is larger than Biomphalaria 
reversa new species and has a sinistral spire with ele¬ 
vated and angled whorls, while B. reversa new species 
has “softer” and rounded whorls. The spiral coefficients 
showed significant differences between both species. 

CONCLUSIONS 

The analysis of the logarithmic spiral of the planispi¬ 
ral shells of the fossil and extant Biomphalaria species 
in this study, combined with shell width and height, 
describes the different spiral growth patterns in different 
species. Also, the analysis is useful as a complementary 
character for species identification, but not significant 
enough to be used by itself. Therefore, to make decisions 
about specimen identification, qualitative moqihologi- 
cal characters had to be used. In fact, the quantitative 
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Figures 15-32. Shells of Biomphalaria species. 15-26. Biomphalaria revena new species, shells in apical, umbilical, and apertural 
view. 15-17. Holotype (FCDPI 7230). 18-20. Paratype a (FCDPI 7830a). 21-23. Paratype b (FCDPI 7830b). 24-26. Paratype c 
(FCDPI 7830c). 27-32. Biomphalaria walteri , shells in apical, umbilical, and apertural view. 27-29. Holotype, Carnegie Museum of 
Natural History, Pittsburgh (CM 103839) (courtesy of Dr. Timothy Pearce). 30-32. FCDPI 6457. Scale bars: 1mm. 


analysis and the traditional qualitative analysis com¬ 
plement each other in the definition of species in the 
genus Biomphalaria. 

A new species is described as a result of the combined 
approaches of statistical analysis and traditional mor¬ 
phological description: Biomphalaria reversa new spe¬ 
cies, from the Queguay Formation (Late Cretaceous) 
of Uruguay. This new species is added to the fresh¬ 
water assemblage of this lithostratigraphic unit, consti¬ 
tuting, with B. walteri, the oldest representatives of 
Biomphalaria worldwide. 
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